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We report measurements of properties of turbulent thermal convection of a fluid with 
a Prandtl number Pr = 4.38 in a cylindrical cell with an aspect ratio T = 0.50. The 
rotational symmetry was broken by a small tilt of the sample axis relative to gravity. 
Measurements of the heat transport (as expressed by the Nussclt number Nu), as well as 
of large-scale-circulation (LSC) properties by means of temperature measurements along 
the sidewall, are present ed. In contradistinc tion to similar experiments using containers 
of aspect ratio T = 1.00 (j Ahlers et aZ.ll200d ) and T = 0.50 (|Chilla et all\2004 ISun et al 



120051: iRoche et aZ.ll2010l) . we see a very small increase of the heat transport for tilt angles 
up to about 0.1 rad. Based on measurements of properties of the LSC we explain this 
increase by a stabilization of the single-roll state (SRS) of the LSC and a de-stabilization 
of the double-roll state (DRS) (it is known from previous work that the SRS has a slightly 
larger heat transport than the DRS). Further, we present quantitative measurements of 
the strength of the LSC, its orientation, and its torsional oscillation as a function of the 
tilt angle. 



1. Introduction 

Turbulent fluid motion driven by a temperature gradient plays an important role in a 
variety of natural and industrial processes. Thus, it has been a topic of i n terest for man y 
decades [for reviews intended fo r a broad audience, see iKadanoffl (l200lh: I Ahlers) (|2009h 
for more specialized reviews, see I Ahlers et all (12009 Lohsc fc Xial (120 



to fu ndamental problems in astro- or geophysics ([Howard &: LaBonte 



(J)]. It is relevant 



1980t iRahmstorl 



20001) . but plays an important role also in applications ra nging from the optimization of 
heat dist ribution insid e aircraft cabins (jKuhn et aZ.ll2009l) to the cooling of giant nuclear 
reactors ( Stacevl[2010l) . 

In order to understand the fundamental processes of turbulent thermal convection, ex- 
periments and simulations focus on simple systems with well defined boundaries known 
as Rayleigh-Benard convection (RBC). Such systems usually consist of a Newtonian fluid 
that is confined by a warm horizontal plate from the bottom and a colder parallel plate 
at a distance L from the top. When the temperature difference AT between the bot- 
tom (T;,) and the top plate (T t ) is not too large, i.e. when the fluid properties do not 
change significantly within that temperat ure range, the Oberbeck-Boussincsq approxi- 
mation (|Oberbecklll879t lBoussinesq||l903l ) can be applied and for a given geometry the 
state of the system is defined by only two dimcnsionless parameters. These are the ratio 
of the driving buoyancy to the viscous and thermal damping forces which is expressed 



f Current address: Department of Physics & Center for the Study of Complex Systems, 
University of Michigan, Ann Arbor, MI 48109 USA 



"Stephan WeissCurrent address: Department of Physics & Center for the Study of Complex Systems, University of Michigc 



by the Rayleigh number 



Ra = 



gaATL 3 



(1.1) 



and the ratio of the kinematic viscosity v to the thermal diffusivity k which is given by 
the Prandtl number 

Pr=-. (1.2) 

K 

Here, g and a denote the gravitational acceleration and the isobaric thermal expansion 
coefficient. 

A major issue of interest in RBC is the heat transport between the bottom and the top 
that is expressed in terms of a dimensionless parameter known as the Nusselt number 



A /; A AT 



(1.3) 



The applied heat current is given by Q, the heat conductivity of the fluid is A, and A is 
the cross sectional area of the cell. Significant progress has been made in understanding 
how he at is transported by the flow and thus, how the Nusselt number depends on Ra 
and Pr (jAhlers et «L||20096| ). Nonetheless there remain unanswered questions. 

In turbulent RBC heat is transported by fluid motion which in part is driven by 
plumes that detach from thermal boundary layers at the bottom or top plate and 
rise or sink due to their bu oyancy (|Kadanofj l200ll IShang et all 120031 ; IXi et all \2004 
Funfschilling fc Ahlerdliool . These plumes are carried by, and by virtue of their buoy- 
ancy in turn drive, a larger coherent flow structure known as the large-scale circulation 
(LSC). For cylindrical cells with aspect ratio r = Z?/L = 1 (D is the diameter) the LSC 
consists of a single convection roll and is known as a single- roll state or SRS. Once in a 
while the SRS collapses but then quickly re-establishes itself. The collapse is known as a 
cessation ( Brown fc Ahlersl 2006 &1 ) . For smaller r (and Pr near five) two counter-rotatini 
rolls, one on top of the other, were found as well (|Xi fc Xiall2008ll Weiss fc Ahlersll2011 
and are known as a double-roll state or DRS . Then the system fluctuates randomly be- 
tween the DRS and the SRS. The fraction of time wsrs spent in the SRS increases 
with increasing Ra, while the fraction wdrs spent in the DRS decreases. It was found 
that the SRS transports heat slightly more efficiently than the DRS, but th e difference 
is only about one to two percent depending on Ra ( Weiss fc Ahlersl |2011&[) . Thus the 
conditional Nusselt number Nusrs of the SRS is larger than Nudrs of the DRS. The 
long-time average, designated simply as Nu, depends on wsrs and wdrs- A bi- modality, 
where for different measurements at the same Ra two different valu es of Nu were foun d 
( Roche et al. 2002 ). or long transients within a single experiment ( Chilla et oi.ll2004h . 



cannot be e xplained by the existenc e of the SRS and the DRS over the parameter range 
explored bv lWeiss fc Ahlers (20116) because they found that the system freely samples 
both states over reasonable time intervals. 

In the past, experiments were done in which the rotational symmetry of the system was 
broken by inclining the cylinder axis at a small angle relative to gravity. Studying these 
slightly anisotropic systems is interesting because many natural convection systems are 
subject to a broken rotational symmetry {e.g. convection on hillsides). But in addition 
one can learn s o meth ing about the LSC and its shape in containers with a vertical axis. 
ICiliberto et al. ( 19961 ) used a sample with a rectangular cross section and Pr ~ 3, and 
tilted their container by 10 degrees in order to fix the orientation of the LSC plane. Within 
their resolution of a percent or so they did not find a change in the heat transport. 
However, they found a significant reduction of the temperature fluctu ations close to 
the bottom plate. Similar experiments were done bv ICioni et al. ( 1997 ) for fluids with 
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Pr ~ 0.025 in a cylindrical sample with V = 1.0. They applied a tilt of 0.5 degrees in order 
to create a preferred orientation of the LSC. They also measured the heat flux through 
the cell for tilt angles up to 4 degrees and found that " [...] the heat transfer is unchanged 
(by no more than 1% even for much higher tilt of 4 degrees), [...]". A more comprehensive 
study of the influence of th e tilt angle on the heat transport in T = 1.0 cylinders was 
done by [AM ers et al\ (|2006[ ) for Pr = 4.4. They found a very small reduction of Nu by 
about 0.5% for a tilt angle of 10°. The changes of the amplitude, and of the dynamics, of 
the LSC were studied as well and found to be substantial. Such a small reduction of the 
heat transport differs from e xperimental r esults reported for cylinders with aspect ratio 
r = 0.5. An investigation bv lChilla et all (|2004l ) found a decrease of Nu of up to 5 % for 
tilt angles as small as /3 = 1.5°. Another set of measurements bv lSun et al. ( 20051 ) also 
found a rather large decrease of Nu, namely by about 2% for a tilt angle of 2°. While these 
two publications report convection with w ater and Pr=4.38, experiments with Helium 
(Pr ~ 1) where done bv lRoche et al. ( 2010h . They made measurements for two different 
tilt angles and found that for tilts of 3.6° the heat transport was reduced by about 1.5% 
in comparison to that for tilt angels of 1.3°. The above results suggest that there is a 
difference between the influence of a tilt on T = 1 cells (almost no reduction of Nu) and 
the effect seen in T = 0.5 experiments (reduction of a few per cent even for very small 
tilt angles). However, our present study for T = 0.50 does not confirm these results, and 
instead finds a very small increase of the heat current. 



Starti ng with the Navier-Stoke s equation and making physically reasonable approxi- 
mations, iBrow^^^hleri (|20086l ) developed a model consisting of two coupled stochastic 
ordinary differential equations that describes the dynamics of the angular orientation and 
the amplitude of the LSC in containers with aspect ratio T = 1. This model was extended 
later to systems with broken ro tational symmetry, including samples tilted relative to 
gravity ( Brown fc Ahlersl 2008a ). Generally good agreement with experimental observa- 
tions was found. Whereas in principle this model should apply to the SRS independent 
of r (albeit with coefficients that depend on V), an extension of the model to the DRS 
would require a larger number of coupled equations and to our knowledge has not yet 
been attempted. 



In this paper, we present a comprehensive experimental investigation of the influence 
of a tilt on the heat transport and the flow structure in a cylindrical RBC sample with 
aspect ratio T = 0.50. In addition to measurements of Nu, we analyze the LSC structure 
and examine the transitions between the SRS and the DRS. One of our findings is that 
the time wsrs spent in the SRS increases with tilt angle while w^rs decreases. This 
result is not surprising because in the SRS gravity can favor both up-flow under the tilted 
top plate and down-flow under the tilted bottom flow, wh ereas for the DRS o nly the flow 
under one of the plates can be favored at the same time ( Chilla et al. 2004f ). For angles 
larger than about 6° the DRS no longer exists and wsrs — 1 ■ We found that the increase 
of wsrs is accompanied by a small increase of Nu. This too is to be expected because 
Nusfjs is large r than Nuns'; (I Weiss fc Ahlersll2011&l). but differs fro m results of previous 
investigations ( Chilla et «dl2004l ; ISun et a/.ll2005l: iRoche et aLll2010l ). W henever possible, 
we com pare our Nu res ults with tho s e of Ahlers et all ( 20061 ) (r = 1.0). lBrown fc Ahlers 
( 2008al) (r = 1.0), and lChilla efall (<2004h (T = 0.5). 



In addition to the Nu measurements, we report on the LSC amplitude, on a Fourier 
decomposition which gives its mode structure, and on the probability distribution of the 
LSC orientation. 
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2. Experimental setup and methods 

We used the medium conve ction apparatus (MCA) described bv lZhong fc Ahlersl (l2010h 



and Weiss fc Ahlersl (20116). Water at an average temperature T m = (T, + T t )/2 



40.00 °C (Pr=4.38) was confined between two copper plates, separated by a distance 
L = 495 mm, and a cylindrical Plexiglas sidcwall of inner diameter D = 247.5 mm, re- 
sulting in an aspect ratio T = 0.500. The sidewall thickness was 6 mm. Three sets, each 
consisting of eight thermistors distributed uniformly in the azinuthal direction, were lo- 
cated at the vertical positio ns z = L/4, L/2 and 3L /4 ab ove the bottom plate. A ll other 
details were as described bv lZhong fc Ahlersl (|2010l ) and IWeiss fc Ahlersl (|20116f ). 

We did experiments for Ra = 1.8 x 10 10 (AT = 3.98 K) and 7.2 x 10 10 (AT = 15.88 K) 
and tilt angles over the range < /3 < 0.12 rad. In a typical run we held AT constant and 
measured the temperatures at all of the thermistors every 3.2 sec for about a day. After 
that, the tilt angle was changed. In addition, for a few measurements we took data for 
very long times (up to 140 hours) to ensure that the statistics does not change for longer 
measurement times. Data that were taken within the first five hours of the experimental 
run were not considered in the data analysis to avoid transient behavior of the system. 

We used the 24 sidew all thermistors to measure the orientation and dynamics of the 



LSC (jBrown et aZ.ll2005i ). To lowest order, the LSC is a flow where warm fluid rises at one 
side of the cell and cold fluid sinks at the opposite. Thus, one can detect the azimuthal 
orientation and the strength of the LSC by fitting a sinusoidal curve of the form 

Tf = T Wi k + 6k cos - 6> fc ^ , (2.1) 

to all eight thermistors at one specific height. Here the index "i" stands for the azimuthal 
location of the thermistors and takes values i = . . . 7. The index "k" denotes the vertical 
location of the thermistor and will in the following have letters "b" (z=L/4), "m" (z=L/2) 
and "t" (z=3L/4). In this way, the amplitude 8k is a measure for the strength of the LSC 
at height level " k" whereas the orientation of the flow is given by the phase 9 k ■ 



3. Results 

3.1. Nusselt number measurements 

Nusselt-number measurements in the absence of a tilt (0 = ) taken with this apparatus 
were published before (see Fig. 3 of Weiss k, Ahlersl ( 201161 )) and agree with previous 



measurements and the predictions of the Grossmann-Lohse model (IGrossmann fc Lohse 



20011) to within a percent or so. For the two Raylcigh numbers investigated in the present 
paper, Nu in the absence of a tilt had the values Nu(0) = 160.4 for Ra = 1.8 x 10 10 
and Nu(0) = 246.4 for Ra = 7.2 x 10 10 . Here we focus on the effect of a tilt and 
normalise the measured Nusselt numbers by Nu(0). In order to improve the accuracy of 
the normalization, we use the reduced Nusselt numbers Nu re d(/?) = Nu(/3)/Ra ' 3 rather 
than Nu. In this way we account for small variations of Ra between different experimental 
runs. 

Figure QJ, shows the normalised Nusselt number Nu re d(/3)/Nu re£ z(0) as a function of 
the tilt angle (3 for Ra = 7.2 x 10 10 . In order to check for hysteretic (history dependent) 
effects, we plot measurement points that were taken while /3 was increased as solid sym- 
bols, whereas measurements with decreasing /3 are shown as open symbols. The data for 
increasing and decreasing tilt agree within their range of scatter and thus no hysteretic 
effect is visible here. 

The behavior of Nu(/3)/Nu(0) differs significantly from measurements reported by 
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Figure 1. The reduced Nusselt number Nu ret i/Nu r ed(0), with Nu red = Nu/Ra 0,300 , as a function 
of the tilt angle f3. (a): Data for Ra = 7.2 x 10 10 and for increasing tilt (solid symbols) and 
decreasing til t (open symbols). The dashed line has a slope equal to the one measured for 
T = 1.00 by lAhlers et al\ (|2006f ). (b): Data for both Ra = 1.8 x 10 10 (squares, red online) 
and Ra = 7.2 x 10 10 (bullets, blue online). Data points in (b) include data for increasing and 
decreasing tilt angles. 



Chilla et all (|2004l) . ISun et all (j2005l) and iRoche et all (|2010l) . We observe a very small 
but well resolved increase of the heat transport with increasing tilt angle. The heat 
transport reaches a maximum near (3 ~ 0.07 rad and decreases for larger [3. Over the 
investigated range of (3 < 0.113 rad Nu is always larger than it is in the horizontal case 
P = 0. However, the increase is very small, with a maximum Nu increase of only about 
0.3%. 

Similar behavior, but with larger scatter of the data, is also observed for Ra = 1.8 x 10 10 
as marked by squares (red online) in figure [Tb. The larger scatter is due to the fact 
that the smaller value of Ra corresponds to a smaller AT where the same temperature 
resolution leads to a larger relative scatter. Due to the scatter of the small-Ra data it 
can not be concluded whether or not there exists a maximum in the range < /3 < 0.12. 
Nevertheless an increase of the heat transport with increasing tilt angles is clearly visible 
in the data. 

The observation of a small increase of heat trans port with increasing tilt seems un- 
expected and is in contrast to previous observation (jChilla et a/.ll2004 ISun et qZ.I|2005 ; 
Roche et a/.ll2010l ). However, we shall show that the reason for this increase can be found 
in a stabilization of the SRS relative to the DRS. Since the SRS transports heat more 
efficiently than the DRS by about one to two percent depending on Ra, a small net 
increase of Nu can be explained. For this purpose we analyze the side- wall measurements 
in the following subsection. 

3.2. Flow-mode transition 



As mentioned in the introduc tion, studies on RBC in cylinders with T = 0.50 (|Xi fe Xia 
20086 ; Weiss fe Ahlersl 20116 ) have shown that in these geometries and for Pr ~ 5, the 
large-scale circulation can exist in form of a single-roll state (SRS) or a double-roll state 
(DRS). The latter consists of two counter-rotating rolls, one on top of the other. The flow 
switches randomly between the two states. The fraction of time that the system spends 
in the DRS depends on Ra, but is in general small in c omparison to the f raction of time 
the system exhibits a S RS. It was shown (see Fig. 19 of IXi fc Xial (|20086l) and Fig. 14 of 
Weiss fc Ahle"r1(|2011&D ) that the heat transport is less efficient when the DRS is present 
and thus the Nusselt number is reduced slightly in comparison to the times when the 
system exhibits the SRS. 

In a slightly inclined convection cell one would expect the SRS to be favored by the 
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Figure 2. (a): The fraction of time wsrs that the system exists in the single-roll state, (b): 
The fraction of time wdr.s that the system exists in the double-roll state, (c): The fraction of 
time wts that the system exists in the transition state. The data are for Ra = 1.8 x 10 10 (red 
squares) and Ra = 7.2 x 10 10 (blue bullets). 



system, since in the DRS the buoyancy adjacent to either the inclined cold (top) or the 
inclined warm (bottom) plat e acts against the direction of the fluid flow, depending on 
the orientation of both rolls ( Chilla et aLll2003) . 

To detect the state of the system, we look at the amplitu des Sk and phases 9t at the 
three different levels z = L/A, L/2 and 3L/4. As described by I Weiss fc Ahlers ( 20116 ). we 
define the system to be in the SRS when \9 t — 0i,\ < 60° and demand that the amplitudes at 
all levels k="t","m" and "b" be larger than 15% of their average values (Sk > 0.15((5fc)). 
On the other hand, we say that the system is in the DRS when \8 t — 6i,\ > 120° and 
5 t ,b > 0.15((5t ! b). When neither the conditions for a SRS nor those for the DRS are 
fulfilled, we call this a transition state (TS). The largest contribution to the TS comes 
from states for which 60 < \0 t - 9 b \ < 120°. 

Figure [2] shows the fraction of time that the system spends in the SRS (a), in the DRS 
(b), and in the TS (c) for the two investigated Rayleigh numbers. The result is clear and 
as expected. With increasing tilt angle, the system spends more time in the SRS and less 
in the DRS or TS. For angles larger than f3 0.06 rad the system consists essentially at 
all times of a single well defined roll, where warm fluid flows along the bottom plate in 
t he uphill directio n and c old fluid flows along the top plate in the downhill direction. 

Weiss fc Ahlers ( 2011&I ) showed that the heat is transported more effectively when the 



system is in the SRS. Therefore, we believe that the preference of the inclined system for 
the SRS explains the increase of Nu with (3. We note, that the location of the maximum 
of Nu(/3) /Nu(0) at f3 0.07 rad roughly coincides with the angle above which the system 
exhibits the SRS during nearly the whole measurement time. Thus, a further increase of 
P cannot result in a further increase of Nu and from there on Nu decreases gradually. To 
illustrate this fact in more detail, we calculate the conditional Nusselt number Nusrs 
that takes only the time intervals into consideration when the system is in the SRS. This 
quantity, divided by Nu, is plotted as a function of (3 in Fig. [3] The plot shows two things. 
First, for no or small tilt angles Nugas is larger than Nu, stating again that the heat is 
transported more efficient in the SRS than in the DRS or the TS. Second, for increasing 
tilt angles Nus rs becomes more nearly equal to Nu since the SRS becomes the dominant 
flow state of the system. 

For cylindrical samples with T = 1.0 the sample is always in the SRS. Consistent 
with this and the phenomenon described above, measurements of Nu did not reveal 
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Figure 3. The conditional Nusselt number Nusks computed only from data taken while the 
system was in the SRS. The results were normalised by Nu and are shown as a function of the 
tilt angle. They are for Ra = 7.2 x 10 10 . The dashed horizontal line marks Nushs/Nu = 1. 



any increase with j3 
corresponding to Nu(/3) 



Instead they showed a very gradual decre ase with increasin g /3, 

~ 20061 ). In 



Nu(0)(l + aj^ u /3) with a^ u 



-0.031 (Ahlers et al 



Fig. Q]a we show a dashed line which has a slope corresponding to the measured aj^ u for 
T = 1.00. The present data for T = 0.50 arc consistent with a similar decease of Nu(/3) in 
the (3 range beyond the maximum where the initial increase due to a diminished presence 
of the DRS no longer occurs. However, the data do not extend to sufficiently large j3 to 
yield quantitative information on this issue. 

In order to compare our results with those of others, we briefly consider l ow-Prandtl- 
number c onvection in cyli n ders wi th T = 0.5. On the basis of measurements bv lAhlers et al 
(20093), IWeiss fc Ahlersl (|20116f ) showed that, for Pr = 0.67, Ra = 1.0 x 10 11 , and 
T = 0.50, the system is in a SRS almost all the time. They found no evidence for a 
DRS. Thus, we expect no increase of Nu with for this case, and by analogy to the 
T = 1.00 measurements a slight decrease of Nu with /? might be expected. Consequently, 
the reduction of the heat transport by about 2% with a chang e of from . 023 to 0.063 
rad in the Ra range from 6 x 10 10 to 2 x 10 11 reported by iRoche et al. ( 201Clh is in 
qualitative agreement with our measurements. However, the reduction of Nu by 2%, cor- 
respo nding to a slope oa j,, ~ —0.5, differs by about a factor of 16 from the measured 



va lue dAhlers et aLll2006D for Pr = 4.3 and T = 1.00. Thus, in order to reco ncile the da ta 
of Roche et al. ( 201dh quantitatively with our measurements and those of Ahlers et 
(l2006h . a very strong Pr dependence of a^y, , would have to be invoked. In view of recent 



numerical calculations ( Bailon-Cuba et al. 20101 E. P. van der Poel et al. 2011 ). which 



revealed that the nature of the LSC has a larger influence on Nu as Pr decreases to or 
below one, such a strong dependence can not be ruled out. 

The reason for the difference be t ween our d ata and th e heat -transport measurements 
as a function of bv lChilla et al. ( 2004 ) and Sun et al. ( 20051 ) (with Pr ~ 5 similar to 
our case) remains unexplained. 



4. Properties of the LSC 

4.1. The average amplitude of the LSC 

Here we investigate the time-averaged strength of the LSC as a function of the tilt angle 
0. Figure shows the dependence of (5k) on for Ra = 7.2 x 10 10 . For no tilt, (5 m ) is 
significantly smaller than (<5 t ) and (5b). This finding is in contrast to results for T = 1.00, 
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Figure 4. The time-averaged amplitude {Sk) of the azimuthal temperature variation along the 
side wall as a function of the tilt angle p. (a): The amplitude for the top {(St), bullets, blue 
online), the middle ((8 m ), diamonds, green online) and the bottom ({5b), squares, red online) 
thermistor row. The data were taken at Ra = 7.2 x 10 10 . (b): The same data as in (a) but 
normalised by their values at P = 0. (c): Comparison of the normalised amplitude at midheight 
(S m (P))/(S(0)) for Ra = 1.8 x 10 1() (squares, red online) and Ra = 7.2 x 10 10 (bullets, blue 
online), (d): Comparison between F = 0.50 data (bullets, blue online, Ra = 7.2 x 10 10 ) and dat a 
for r = 1.0 [up-pointing triangles, yellow online, Ra = 2.8 x 10 9 , fro mlBrown fc Ahlersl (12008 of ); 
and down-pointing triangles, purple online, Ra = 9.4 x 10 10 , from lAhlers et al.\ (|2006f )]. The 
solid lines in (a) and (d) are fits of Eg. 14. II to the data. The coefficients derived from those fits 
are shown in Table |4~T1 



where (6 m ) was larger ( Brown fc Ahlersl 120071) . but in agreement with previous results 
for T = 0.50 dWeiss fc Ahlersl 1201161) . With increasing tilt angles all three (4) increase. 
Such an increase was also found for T = 1, where however only the temperature a t 
the horizontal mid-plane ((6 m )) was investigated as a function of /? ( Ahlers et a^.ll2006h . 
The increase of the (4) is smooth and does not show any obvious change at or beyond 
f3 ps 0.06 where the DRS and TS cease to exist. 

The increase of (4) with /3 is not linear. The slope d (5k)/d(3 decreases w ith increasing 



p. This is also in agre ement with the re s ults for T = 1 ( Ahlers et al. 20061 ). Based on the 



model for the LSC bv lBrown fc Ahlersl (2008o|), we fit the sinusoidal function 



(4) = (4(0)>-[l + asin(/3)] (4.1) 

to the data for a quantitative comparison. The corresponding fits are shown as solid lines 
in figure H^,. The coefficients (4(0)) and a are listed in table 14. ll The coefficient a does 
not depend significantly on the level. The amplitudes, normalised by their values (4(0)) 
for /3 = 0, are plotted in figure|4jD. One sees that the relative increase of (4) is the same 
for all levels. We shall show below in Sec. 14.41 that this is the case also for the probability 
distributions of the (4)- 

To compare the strength of the LSC for both Ra values that we investigated, we plot 
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Figure 5. The ratio of the energy in the first Fourier mode Ei. m to the total energy of all four 
modes E to t for Ra = 1.8 x 10 10 (squares, red online) and for Ra = 7.2 x 10 10 (bullets, blue 
online). The measurements here are for the middle thermistor row. Measurements for the top 
and the bottom give nearly the same result. 



in Fig. dfc the amplitude at mid-height (<5 m (/3))/(<5 m (0)} for Ra = 1.8 x 10 10 (squares, red 
online) and for Ra = 7.2 x 10 10 (bullets, blue online). One sees that for the larger Ra the 
curve starts to bend earlier. This might suggest that a simple sinusoidal fit is not sufficient 
to characterise (Sk((3)), and that higher-order terms might be needed. Nonetheless, for 
the present purpose we shall retain the use of Eq. 14.11 

In Fig. SJi and Table 14.11 we compare our m easurements of (S m ) (T = 0.50) with 
measureme nts for T = 1.00 a t Ra = 2.8 x 10 9 ( Brown fc Ahlers! 12008 al ) and at Ra = 
9.4 x 10 10 (jAhlers et a/.ll2006l ). One sees that at nearly the same Ra the increase of (S m ) 
is stronger (the coefficient a is larger) for V = 0.50 than it is for V = 1.00. The data also 
suggest a significant decrease of a as Ra increases at constant T. 

4.2. Higher-order modes 

In the previous section (Sk) was obtained from a fit of Eq. l2.ll to the side- wall temperature- 
measurements, and thus it represents the amplitude of the first Fourier mode of the 
azimuthal temperature variation. To gain information about the shape of the LSC it is 
useful a lso to compare the e nergy E\_k = Si of this lowest mode with that of the higher 



modes (|Stevens et all 1201 it IWeiss fc Ahlersl 12011 at ) . Since the behavior is quite similar 
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Figure 6. The time-averaged frequencies at which the amplitudes 5k drop below 15% of their 
average values, (a): Frequency of events for the top (bullets, blue online), the middle (diamonds, 
green online) and the bottom (squares, red online) thermistor row. Measurements were taken 
at Ra = 7.2 x 10 10 . (b): Frequency of events for the middle thermistor row for Ra = 1.8 x 10 10 
(triangles, red online) and Ra = 7.2 x 10 10 (stars, blue online). 



for all three thermistor rows, we discuss only the analysis for the middle row. Figure [5] 
shows the ratio between the energy £a >TO of the first mode to the total energy E to t,m of 
all four accessible modes (we have eight thermistors at one height and thus have access 
to four Fourier modes). The data show an increase of Ei tTn /E tot m with increasing tilt 
angle. One sees that higher-order modes are suppressed by increasing the tilt. We also see 
that Ei^m/ E to t,m is smaller for Ra = 1.8 x 10 10 than it is for Ra = 7.2 x 10 10 . However, 
as P increases, E\_ m j E to t,m reaches similar values, larger than 0.9, for both Ra. 

4.3. Events 

An event is defined as a drop of the amplitude 6k below 15% of its time-averaged value 
((5fc). In this section we report on the time-averaged event frequency cuf. as a function of 
the tilt angle. The cause of events may vary and depends on T. While for T = 1.0 most of 
the events are cessations [almost no flow-mode transitions are found dXi fc Xiall2f)08feh ]. 



most of the events for T = .50 are caused by flow-mode transitions [cessations are very 



rare ([Weiss fc Ahlersll2011M) ]. 

Figure [5^ shows the event rate for the three different heights as a function of the tilt 
angle (3. The data points show considerable scatter because the number of events that 
occurred during the measurement times of about a day for ea ch point was small. Fo r 
j3 = the event rate is comparable to previous measurements ( Weiss fc AhlerslboilSl ). 



No significant difference can be observed between the event rate of the top, the middle, 
and the bottom thermistor row. This can be understood since these events are mostly 
flow-mode transitions at which a smaller roll appears at the top or the bottom, grows in 
size and replaces the original roll. During such transitions a dip (an event) can be observed 
in all of the three amplitudes (note however, that this is also true for cessations). 

The decrease of uj 1 ^ with increasing tilt angle, and the observation that for /3 > 0.06 
no events occurred, is consistent with the results for the fraction of time wsrs that 
the system spent in the SRS as reported above in Fig. [5^. As the SRS becomes more 
dominant and the DRS becomes rare, flow-mode transitions (i.e. transitions between 
these two states) become less frequent. When, for j3 > 0.06 only the SRS is found, 
no transitions can occur. The finding of a reduction of event s with increasing tilt is in 



accordance to the results for V — 1.00 ([Brown fc Ahlersll2008al ). There however, the event 
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Figure 7. The probability-density functions of the normalised LSC strength 8k /{5k) for the top 
(bullets, blue online), the middle (diamonds, green online) and the bottom (squares, red online) 
thermistor row. Shown are data for (a) /? = and (b) j3 = 0.035 rad. The solid lines (green 
online) are Gaussian fits to the right side of the peak of S m /{5 m )- The experiments were done 
at Ra= 7.2 x 10 10 . 



rate was already significantly smaller in the horizontal case (1.7 per day) and was mainly 
caused by cessations. 

In Fig.|BjD we compare the event rates of the middle thermistor row uj^ for experiments 
with Ra = 1.8 x 10 10 and Ra = 7.2 x 10 10 . For j3 < 0.05 the data suggest a slightly smaller 
event rate for the smaller Ra, but due to the scatter a definitive conclusion can not be 
drawn. 

4.4. The probability distribution of Sk 

We showed previously ( Weiss fe Ahlersl 1201161 ) that the probability-density functions 
(PDFs) p(Sk) of the amplitudes 4 at different levels k collapse when the 4 are scaled by 
their corresponding time averaged values {Sk)- This is well illustrated by figure [7£i where 
for the horizontal case (j3 = 0) the PDF's of the normalised amplitudes are plotted. 
The graph looks very similar to the one shown bv IWeiss fe Ahlersl ( 2011 feh . The PDF's 
are asymmetric. At 5k/{Sk) = they start with a finite slope that increases until an 
inflection point is reached. From there on p{5k/ {Sk)) continues to grow until it reaches a 
maximum. The maximum and the right tail of the PDF can be described fairly well by a 
Gaussian distribution as indicated by the solid lines in figurc[7l All this is consistent with a 
description of Sk as a st ochastically drive n amplit ude diffusing in an asymmet ric potential 
which was derived by Brown fc Ahlersl ( 20086 ) and extended recently by Assaf et al. 
(|201ll ). albeit with parameters for T = 1.00. 

Figure [TJd shows data similar to those in (a) , but here the cell was tilted by j3 = 
0.035 rad. The effect on p(Sk/{8k)) is clearly visible. The peak becomes significantly 
narrower (and thus higher) than it is in the horizontal case. Since also the slope at 
Sk/{Sk) = becomes smaller, the peak becomes more nearly symmetric. This result is 
in accordance with the reduction of flow-mode transitions and the stabilization of the 
SRS since flow- mode transitions and other events (see section I4.3|) result in a temporary 
decrease of Sk, and thus in a skewed distribution. 

For a quantitative analysis, we calculate the standard deviation 



a k = ((4/(4) - i) 2 ) 1/2 



(4.2) 



of the normalised amplitude for the middle thermistor row k = m. We note that a similar 
calculation for the other thermistor rows gives very similar results, as indicated by figure 
[7] In Fig. [8]we show the inverse 1 jo m as a measure of the stability and the coherence of the 
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Figure 8. Inverse of the standard deviation a m of the amplitude distribution as a function of the 
tilt angle /3. (a): Comparison between Ra = 1.8 x 10 10 (squares, red online) and Ra = 7.2 x 10 10 
(bullets, blue online), (b): Comparison between F = 0.50 (b ullets, blue online , Ra = 7 .2 x 10 10 ) 
and T = 1.00 (triangles, yellow online, Ra = 2.8 x 10 s , from lBrown fc Ahlersl l |2008aft ). 



LSC. It is interesting that in the horizontal case l/er m is the same for both investigated 
Ra values. However, when the sample is tilted, l/o~ m increases somewhat faster for the 
larger Ra. 

In Fi g. |8Jd we compare our data (r = 0.50) with those for T = 1.00 of iBrown fc Ahlersl 



( 2008q| ). Note, that there is a difference of almost three orders of magnitude in Ra. It 
is not a surprise to see that in the horizontal case the flow is more stable (larger 1/er) 
for r = 1.00 than it is for T = 0.50. This was already found and discussed in previous 



publications ( Xi fc Xia 2008fcl fal; IWeiss fc Ahlersl 20116 ). However, it is surprising that 



for p > 0.8 the T = 1.00 data show the same slope as the T = 0.50 data did in the 
/3 range that we investigated. This seems remarkable, especially since the Ra values at 
which the two data sets were taken differ by two and a half orders of magnitude. 

4.5. The orientation of the LSC 
By inclining the cell from its horizontal position, the rotational symmetry is broken and 
the azimuthal orientation 6k of the LSC plane is no longer randomly distributed but 
tends to align in the direction of the tilt. This process can be investigated by looking at 
the probability-density function of 6 k at the different heights. 

4.5.1. The horizontal case f3 = 

For the perfectly horizontal and rotationally invariant system all angles 8k should be 
sampled equally and the (normalised) PDF p(6k), < 6k < 2tt, should be a constant 
equal to l/(27r). PDFs for 6k at Ra = 7.2 x 10 10 are shown in Fig.[!|]for (a) the horizontal 
case (3 = 0.000 and (b) a tilt of (3 = 0.035 rad. One sees that already for the horizontal case 
there is a small anisotropy. The circulation plane is less likely to be found near 6/2tt w 0.3 
and more likely to be found near 6/2ir —0.2. The PDF has to be 27r-periodic. It turns 
out that the simplest function with this property, namely 

p(6k) = 1/(2tt) + A k ■ sm(6 k - fa) , (4.3) 

fits the data very well. Here the phase fa, although it should be equal at all k, is arbitrary 
since it depends on the arbitrarily selected origin of the azimuthal coordinate system 
which may vary from one experimental setup to another. A fit of Eq. 14.31 to the data for 
k = t (the top thermistor row, blue bullets) is shown in Fig.[5£i as a solid line (blue online). 
The fits for k = m (solid diamonds, green online) and k = b (solid squares, red online) 
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Figure 9. Probability density functions of the orientations 6^ of the LSC plane for 
Ra = 7.2 x 10 10 at tilt angles of (a) /3 = 0.000 and (b) /3 = 0.035 rad. The data are for 
the top (bullets, blue online), the middle (diamonds, green online) and the bottom (squares, red 
online) thermistor row. The open gr een circles in (a) show p(0 m ) from a previous experiment 
with Ra = 9.0 x 10 10 and /3 = 0.000 (|Weiss fc Ahlerslboilftft . The lines in (a) are fits of Eg. [4731 
to the data, and the lines in (b) are fits of Eg. 14. 41 to the data. 



are not shown for clarity, but they are equally good. The amplitude Ak is a measure of 
the extent of the deviations from the expected uniform distribution of the rotationally 
invariant system. For the example of Fig. |H] we found Ab = 0.068, A m = 0.056, and 
At = 0.057, all with probable errors of 0.003. As expected, there is little if any significant 
height dependence. 

Also shown in Fig. [SJi, as op en circles (green onlin e), are results from a previous 
investigation for Ra = 9.0 x 10 10 (jWeiss fc Ahlerdl20116l) ). The dashed line ( green online) 
in Fig. |9K is the fit of Eq. 14.31 to those data (in this case the phase <j) m is different because 
the arbitrary origin of the azimuthal coordinate was different). That fit gave A m = 0.044, 
which is somewhat lower than but close to the present resu lt. A third data set, for T = 0.50 
and Ra = 5.7 x 10 10 , was reported by IXi fc Xial (|20086h and suggests similar values of 
the Ak- 

The inhomogencity of p(0k) might be due, for instance, to small anisotropics in the 
experiment, such as for example of the spatial temperature di stribution in the top an d 
bottom plates, or to small inhomogcneitics of the sidcwall (jBrown fc Ahlersl 12008 al ). 
However, in experiments with T = 1.0 the anisotropy in p(6k) is significantly stronger 
and even for the horizontal case the PDF consists of a well defined peak w hich can be fit 
well by a Gaussian distribution ( Brown fc Ahlersl 2006a : IXi fc Xial l 2008 q| ). It was shown 
that the Coriolis force caused by the Earth's rotation leads to a preferred orientation of 



the LSC (IBrown fc Ahlersl 2006a ) with a PDF in agreement with the measurement. As 
pointed out by Xi fc Xia ( 20080 ). according to the model of Brown fc Ahlersl ( 2006a ) 
the influence of the Coriolis force should be even stronger for V = 0.50 than it is for 
T = 1.00 because the vertical part of the LSC, which causes the preferred orientation, 
is longer. However, the dynamics of the LSC for T = 0.50 is more erratic and a SRS is 
often destroyed by events as documented above in Sec. I4.3l and Fig. [6] (e. g., by flow mode 
transitions) and r eappears within short time intervals, often at a different orientations 
( Xi fc Xia 2008 al ). This tends to randomize the orientation, counter-acts the orienting 
influence of the Coriolis force, and leads to a broader PDF of 8k ■ 
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Figure 10. The inverse of the standard deviation ag of the distributions as a function of the 
tilt angle ft. (a): Comparison between Ra = 1.8 x 10 10 (squares, red online) and Ra = 7.2 x 10 10 
(bullets, blue online). The solid lines are fits of the straight line l/cr$ — p + q- ft to the data. The 
fits yielded the offsets and slopes p = 1.0 ± 0.1, q = 67 ± 2 for Ra = 1.8 x 10 10 and p = 1.1 ± 0.1, 
q = 61 ± 2 for Ra = 7.2 x 10 10 . The insert shows the same data, but with ag (rather than its 
inverse) plotted against ft. (b): The bullets (blue online) are fo r Ra = 7.2 x 1 l0 . T he triangles 
(yellow online) are data for V = 1.00 at Ra = 9.43 x 10 io from lAhlers et al\ ([20061') (these data 
were multiplied by 2tt in order to correct a mistake in Fig. 12 of lAhlers et al. "( |2006l 'l'l. 



4.5.2. The inclined case j3 > 

The difference between the horizontal case (Fig. |9K) and the inclined one (Fig. [9L) 
is clearly visible. At ft = 0.035 rad only orientations close to 6/2n = 0.13 existed. This 
shows that the orientation of the LSC is concentrated near one azimuthal position. It 
also suggests that the single roll is more persistent and is less often destroyed by events 
(see Fig.[6|), such as flow- mode transitions or cessations, which would tend to randomize 
the orientations. 

To quantify the narrowing of the PDF, we fit the Gaussian distribution 

p{e m ) = Ae^(- {6m -^ o) ^ (4.4) 

to p{9 m ). Although this function is not 2-7T periodic (as in principle p{0k) should be), 
it provides an excellent approximation when the peak is narrow and the tails do not 
extend too far. Nonetheless, for the fit we only used points close to the maximum since 
for small tilt angles one tail of the Gaussian is usually superimposed upon the other 
due to the 27T-periodicity of 6 m . The insert in figure [TDk shows the decrease of ag with 
increasing tilt angle. However, by plotting its inverse 1/ag in figure [TU] one sees that this 
quantity follows a linear trend in the range of tilt we have investigated. This trend seems 
to depend only slightly on Ra. A fit of the straight line 1 jag = p + q ■ (3 to the data 
revealed a slightly larger slope for the smaller Ra (q = 67 ± 2 for Ra = 1.8 x 10 10 and 
q = 61 ± 2 for Ra = 7.2 x 10 10 ). 

4.6. The oscillation of the LSC and its relation to the Reynolds number 

In addition to the rand om diffusion of the orientation and the amplitude of the LSC 
( Brown fc Ahlersl 12008 bh. the LSC also performs a torsional oscillation when it is in the 



SRS [Funfschilling fc Ahlersl (|2004l )]. The orientations of a single roll oscillate with the 
same frequency near the top and the bottom of the roll, but with a phase shift of it relative 
to each other. This motion is especially interesting since, for T = 1, the oscillation period 
is the same as the turnover time of the single roll. Thus it can be used to calculate 
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the Re ynolds number Re of the flow (|Qiu fe Tongll2002t iBrown et alll2007t lAhlers etal 
2009 fef ). 



Torsional oscillati ons were studied experimentally for T = 1 ([Brown et "Zl l2007|: Funfschilling et al 
20081: |Xi et al.ll2009l) . They could be explained in terms of a stochastic model bv I Brown fe Ahlersl 
(|2009i r For r = 0.5 there are stronger amplitude fluctuations and disruptions by flow- 
mode transitions than for T — 1.00 and thus these oscillations are harder to detect 
experimentally. Nevertheless, an analysis using only time segments for which the system 
was in the S RS and averaging over m any of these revealed that also for this T a torsional 



mode exists (Weiss fc Allien 



201161). 



As w as done before ( Funfschilling fc Ahlers 2004 ; Ahlers et al. 2006 ; Weiss fc Ahlersl 



2011a ). we uses the correlation functions 



ki ,&2 



(4.5) 



of the differences between the orientations of the LSC at different heights in order to 
detect the torsional oscillations of the single roll. Here (•) indicates a time average, and 



ovk^kj = o kl - a 

where each of the ki, kj (i,j = 1, . . . , 4) can be "6" 
Eq. 14.51 is then normalised according to 



C 



ki ,k 2 



(4.6) 

or "£" . The function given by 

(4.7) 



We focus here on the auto-correlation function of the difference 89t,b and the cross- 
cor relation function betwe en the differences 89 t _ m and 89b^ m . As explained in detail 
by Weiss fc Ahlersl ( 20116 ). we calculated the correlation functions only for the time- 
segments during which the system was in the SRS and then averaged over all such 
segments. Figure [TTI shows the auto- and the cross-correlation functions {C\' b and C t b '™) 
for different tilt angles. For the horizontal case, one can detect a clear anti-correlation 



for r = in C°'™ and an oscillating behavior that manifests itself both in Cf'™ and C\ b . 
However, both correlation functions decay rather quickly, which is due to the erratic 
nature of the LSC associated with the randomly occurring flow-mode transitions. 

With increasing tilt angle, the twist oscillations of the correlation functions become 
more obvious because the oscillation period becomes shorter and the decay time becomes 
longer. In addition, the correlation functions become smoother because the lengths of the 
SRS-segments in the time series increase with increasing tilt angle (see Sec. 1 .3 . 2 [) . 

It is desirable to quantify the change of the torsional oscillation-frequency with increas- 
ing j3. A fit of an appropriate oscillating but decaying function to the data turned out 
to be difficult because of the rapidly decaying envelope and the experimental irregulari- 
ties of the functions that are caused by the relatively short time periods over which the 
SRS was available, particularly at the smaller (i. Thus, we used two ways to determine 
the oscillation period T p . We measure the distance between the two maxima closest to 
t = in the cross-correlation function C t '™, and we determine the distance between 

the first two minima in the auto-correlation function C\' b . Both values are not exactly 
the same because of the decaying envelope which shifts the maxima (minima) slightly 
towards smaller values of r. The results obtained this way are shown in figure IT2k. One 
sees that the periods calculated from C\\ (bullets, blue online) are slightly larger than 

those determined from C t b '™ (squares, red online). However, both follow the same trend, 
namely a decrease of the period with increasing tilt. 



-,/./< 
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Figure 11. The auto-correlation function C\\ (circles, blue online) and the cross-correlation 
function C t '™ (squares, red online) as a function of the delay time t. Shown are correlation 
functions for different tilt angles /3. All experiments were conducted at Ra = 7.2 x 10 10 . 
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Figure 12. (a): Torsional periods as measured from C\'™ (open squares, red online) and Cj'£ 

(open circles, blue online) for Ra = 7.2 x 10 10 . (b): Averages T p of the torsional periods obtained 
by the two methods for Ra = 1.8 x 10 10 (solid squares, blue online) and Ra = 7.2 x 10 10 (solid 
circles, blue online), (c): Normalised torsional frequency w(/3)/w(0) for Ra = 1.8 x 10 10 (solid 
squares, red online) and for Ra = 7.2 x 10 10 (solid circles, blue online). 

In Fig. [T2"b we display the averages T p of the torsional periods obtained by the two 
methods for Ra = 1.8 x 10 10 (solid squares) and Ra = 7.2 x 10 10 (solid circles) as a function 
of /3. They can be represented quite well by the straight lines T p = r + s ■ (3, as shown by 
the lines in the figure. Least-squares fits yielded r = 143 (215) s and s = —419 (—511) 
s/rad for Ra = 7.2 x 10 10 (Ra = 1.8 x 10 10 ). 

Figure [12b shows the normalised oscillation frequencies U)(f3)/u>(0) = T(0)/T(/3). We 
normalised the results in order to compare the increase of the frequency for Ra = 1.8 x 
10 10 with the one for Ra = 7.2 x 10 10 . The /3-dependence of the normalised frequency 
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is similar for both Ra; however, there is a slightly smaller increase for the smaller Ra. 
If we assume that the oscillation f requency is propo rtional to the average velocity of 



the flow (i.e. the Reynolds number (IQiu fc Tonal2002l) k then the results here should be 



comparable with the increase of the LSC amplitude Sk as discussed in Sec. 14.11 Although 
there is a qualitative similarity, one sees that 5(f3)/5(0) increases more rapidly with j3 
than does ui((3)/ui(Q), by about a factor of three or so. Another difference is that for 
w(/3)/w(0) the data for the smaller Ra show a slightly smaller increase in comparison 
with the large-Ra data; for S((3)/5(0) it is the other way around. The reasons for these 
differences are unclear. On the one hand the mapping between the azimuthal temperature 
variation that yields S and the velocity field may not be as simple as assumed. On the 
other hand, the LSC oscillation frequency for T = 0.50 may not simply be proportional 
to the Reynolds number of the SRS. 



5. Conclusion 

In this paper we reported on an investigation of the effect of a small tilt on the heat 
transport and the flow structure in turbulent Raylcigh-Bcnard convection (RBC) in a 
cylinder with aspect ratio V = 0.50. In particular we showed results for tilt angles up to 
j3 = 0.122 rad and for the two Raylcigh numbers Ra = 1.8 x 10 10 and Ra = 7.2 x 10 10 . 

We found that a small tilt enhanced the heat transport very slightly. This differs from 
the case of RBC in cylindrical samples with aspect ratio T = 1.0 where a very small 
decrease had been observed. By looking at the ratio of the times that the system spends 
in the single- rolls state (SRS) and the double- roll state (DRS), we found that the SRS 
is stabilized by the tilt and the time that the system spends in the DRS is reduced. 
For /3 > 0.06 the SRS is found essentially all the time. Since the SRS transports heat 
more effectively than the DRS, its stabilization explains the slight increase of the Nusselt 
number Nu with increasing tilt. 

To investigate the behavior of the large-scale circulation (LSC) in the presence of a tilt, 
we measured the average amplitude of its first Fourier mode (dk) and found an increase 
with increasing tilt. The relative increase was slightly larger for smaller Ra. A stabilisation 
of the flow could be seen as well in the event rate, which is the frequency at which the 
amplitude drops below a critical value taken to be 15% of its average value. Decreasing 
numbers of events were found with increasing tilt. For tilt angles larger 0.06 rad no events 
could be detected during the experimental run time of a day for each data point. 

Also the dynamics of the LSC orientation 8k was influenced by the tilt. The orientation 
was concentrated in one direction and random diffusion over all angles was suppressed. 
This could be seen in the probability density function of 0k which became narrower. 
Over the investigated range of tilt angles, the inverse of the standard deviation (1/crg) of 
p(9k) was proportional to the tilt angle. The torsional oscillation at which the top part 
of the single roll oscillates relative to the bottom part became more pronounced and the 
oscillation frequency increased with increasing tilt angle. For T = 1.00 the frequency had 
been shown to be a measure for the Reynolds number. Thus its increase indicates an 
increasing flow velocity. 

The stabilisation of the LSC can also be seen in a Fourier analysis of the first 4 modes 
that can be measured from the azimuthal temperature profile along the sidewalk With 
increasing tilt the relative energy in the first mode increases at the expense of the others. 
While in the horizontal case the relative energy in the first mode is smaller for smaller 
Ra, for increasing tilt the energy of the first mode becomes nearly independent of Ra. 



IBHephan WeissCurrent address: Department of Physics & Center for the Study of Complex Systems, University of Michigc 

This work was supported by the U.S. National Science Foundation through Grant 
DMR-1158514. 



REFERENCES 
Ahlers, G. 2009 Turbulent convection. Physics 2, 74-1-7. 

Ahlers, G., Bodenschatz, E., Funfschilling, D. & Hogg, J. 2009a Turbulent Rayleigh- 
Benard convection for a Prandtl number of 0.67. J. Fluid Mech. 641, 157-167. 

Ahlers, G., Brown, E. & Nikolaenko, A. 2006 The search for slow transients, and the effect 
of imperfect vertical alignment, in turbulent Rayleigh-Benard convection. J. Fluid Mech. 
557, 347-367. 

Ahlers, G., Grossmann, S. & Lohse, D. 20096 Heat transfer and large scale dynamics in 
turbulent Rayleigh-Benard convection. Rev. Mod. Phys. 81, 503-538. 

Assaf, M., Angheluta, L. & Goldenfeld, N. 2011 Rare fluctuations and large-scale circu- 
lation cessations in turbulent convection. Phys. Rev. Lett. 107, 044502. 

Bailon-Cuba, J., Emran, M. & Schumacher, J. 2010 Aspect ratio dependence of heat 
transfer and large-scale flow in turbulent convection. J. Fluid Mech. 655, 152-173. 

BOUSSINESQ, J. 1903 Theorie analytique de la chaleur, Vol. 2. Paris: Gauthier-Villars. 

Brown, E. & Ahlers, G. 2006a Effect of the Earth's Coriolis force on turbulent Rayleigh- 
Benard convection in the laboratory. Phys. Fluids 18, 125108-1-15. 

Brown, E. & Ahlers, G. 20066 Rotations and cessations of the large-scale circulation in 
turbulent Rayleigh-Benard convection. J. Fluid Mech. 568, 351-386. 

Brown, E. & Ahlers, G. 2007 Temperature gradients, and search for non-Boussinesq effects, 
in the interior of turbulent Rayleigh-Benard convection. Europhys. Lett. 80, 14001-1-6. 

Brown, E. & Ahlers, G. 2008a Azimuthal asymmetries of the large-scale circulation in tur- 
bulent Rayleigh-Benard convection. Phys. Fluids 20, 105105-1-15. 

Brown, E. & Ahlers, G. 20086 A model of diffusion in a potential well for the dynamics 
of the large-scale circulation in turbulent Rayleigh-Benard convection. Phys. Fluids 20, 
075101-1-16. 

Brown, E. & Ahlers, G. 2009 The origin of oscillations of the large-scale circulation of tur- 
bulent Rayleigh-Benard convection. J. Fluid Mech. 638, 383-400. 
Brown, E., Funfschilling, D. & Ahlers, G. 2007 Anomalous Reynolds-number scaling in 

turbulent Rayleigh-Benard convection. J. Stat. Mech. 2007, P10005-1-22. 
Brown, E., Nikolaenko, A. & Ahlers, G. 2005 Reorientation of the large-scale circulation 

in turbulent Rayleigh-Benard convection. Phys. Rev. Lett. 95, 084503-1-4. 
Chilla, F., Rastello, M., Chaumat, S. & Castaing, B. 2004 Long relaxation times and 

tilt sensitivity in Rayleigh-Benard turbulence. Euro. Phys. J. B 40, 223-227. 
Ciliberto, S., ClONl, S. & Laroche, C. 1996 Large-scale flow properties of turbulent thermal 

convection. Phys. Rev. E 54, R5901-R5904. 
Cioni, S., Ciliberto, S. & Sommeria, J. 1997 Strongly turbulent Rayleigh-Benard convection 

in mercury: comparison with results at moderate Prandtl number. J. Fluid Mech. 335, 111— 

140. 

E. P. van der Poel, Stevens, R. & Lohse, D. 2011 Connecting flow structures and heat 

flux in turbulent Rayleigh-Benard convection. Phys. Rev. E 84, 045303(R)-1 - 4. 
Funfschilling, D. & Ahlers, G. 2004 Plume motion and large scale circulation in a cylindrical 

Rayleigh-Benard cell. Phys. Rev. Lett. 92, 194502-1-4. 
Funfschilling, D., Brown, E. & Ahlers, G. 2008 Torsional oscillations of the large-scale 

circulation in turbulent Rayleigh-Benard convection. J. Fluid Mech. 607, 119-139. 
Grossmann, S. & Lohse, D. 2001 Thermal convection for large Prandtl number. Phys. Rev. 

Lett. 86, 3316-3319. 

Howard, R. & LaBonte, B. 1980 The sun is observed to be a torsional oscillator with a period 

of 11 years. Astrophys. J. 239, L33-L36. 
Kadanoff, L. P. 2001 Turbulent heat flow: Structures and scaling. Phys. Today 54 (8), 34-39. 
Kuhn, M., BOSBACH, J. & Wagner, C. 2009 Experimental parametric study of forced and 

mixed convection in a passenger aircraft cabin mock-up. Building and Environment 44 (5), 

961 - 970. 



Turbulent RBC in a tilted T = 0.50 cylinder 



19 



Lohse, D. & XlA, K.-Q. 2010 Small-scale properties of turbulent Rayleigh-Benard convection. 

Annu. Rev. Fluid Mech. 42, 335-364. 
Oberbeck, A. 1879 Uber die Waxmeleitung der Fliissigkeiten bei Beriicksichtigung der 

Stromungen infolge von Temperaturdifferenzen. Ann. Phys. Chem. 7, 271-292. 
QlU, X. L. & Tong, P. 2002 Temperature oscillations in turbulent Rayleigh-Benard convection. 

Phys. Rev. E 66, 026308-1-12. 
Rahmstorf, S. 2000 The thermohaline ocean circulation: A system with dangerous thresholds? 

Climate Change 46, 247-256. 
Roche, P. E., Castaing, B., Chabaud, B. & Hebral, B. 2002 Prandtl and Rayleigh numbers 

dependences in Rayleigh-Benard convection. Europhys. Lett. 58, 693-698. 
Roche, P.-E., Gauthier, F., Kaiser, R. & Salort, J. 2010 On the triggering of the ultimate 

regime of convection. New J. Phys. 12, 085014. 
Shang, X. D., Qiu, X. L., Tong, P. & Xia, K.-Q. 2003 Measured local heat transport in 

turbulent Rayleigh-Benard convection. Phys. Rev. Lett. 90, 074501. 
Stacey, W. M. 2010 Fusion: An Introduction to the Physics and Technology of Magnetic Con- 
finement Fusion. New York: Wiley. 
Stevens, R.J. A.M., Clercx, H.J.H. & Lohse, D. 2011 Effect of plumes on measuring the 

large scale circulation in turbulent Rayleigh-Benard convection. Phys. Fluids 23, 095110- 

1-11. 

Sun, C, Xi, H. D. & XlA, K. Q. 2005 Azimuthal symmetry, flow dynamics, and heat transport 
in turbulent thermal convection in a cylinder with an aspect ratio of 0.5. Phys. Rev. Lett. 
95, 074502-1-4. 

Weiss, S. & Ahlers, G. 2011a The large-scale flow structure in turbulent rotating Rayleigh- 
Benard convection. J. Fluid Mech. DOI:10.1017/jfm. 2011. 392, pub. online Nov. 1 2011. 

Weiss, S. & Ahlers, G. 20116 Turbulent Rayleigh-Benard convection in a cylindrical container 
with aspect ration T = 0.50 and Prandtl number Pr=4.38. J. Fluid Mech. 676, 5-40. 

Xi, H. D., Lam, S. & XlA, K.-Q. 2004 From laminar plumes to organized flows: the onset of 
large-scale circulation in turbulent thermal convection. J. Fluid Mech. 503, 47-56. 

Xi, H.-D. & XlA, K.-Q. 2008a Azimuthal motion, reorientation, cessation, and reversals of the 
large-scale circulation in turbulent thermal convection: A comparative study in aspect ratio 
oner and one-half geometries. Phys. Rev. E 78, 036326-1-11. 

Xi, H.-D. & XlA, K.-Q. 20086 Flow mode transition in turbulent thermal convection. Phys. 
Fluids 20, 055104-1-15. 

Xi, H.-D., Zhou, S.-Q., Zhou, Q., Chan, T.-S. & Xia, K.-Q. 2009 Origin of the temperature 
oscillation in turbulent thermal convection. Phys. Rev. Lett. 102, 044503-1-4. 

Zhong, J.-Q. & Ahlers, G. 2010 Heat transport and the large-scale circulation in rotating 
turbulent Rayleigh-Benard convection. J. Fluid Mech. 665, 300-333. 



